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ABSTRACT 

An  IBM  7090  computer  program  is  described  which  calculates  the  distri¬ 
bution  of  distances  between  the  point  of  activation  of  a  weapon  and  a  moving 
target  submarine  taking  into  account  the  estimated  component  attack  errors. 
The  model  used  is  more  flexible  and  realistic  than  similar  past  efforts  and  is 
expected  to  produce  more  reliable  submarine  kill  probabilities.  The  miss 
distances  are  computed  by  Monte  Carlo  simulation  of  the  actual  tracking  and 
firing  tactics.  They  are  plotted  by  an  SC  4020  plotter,  first  in  ascending 
order,  then  as  a  cumulative  frequency  distribution.  Flow  charts,  a  listing  of 
the  FORTRAN  program,  and  a  sample  calculation  are  included. 
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I.  INTRODUCTION 


The  determination  of  the  probability  of  killing  un  evading  submarine  has 
long  been  a  major  problem  of  antisubmarine  warfare.  Past  attack  models  have 
been  usually  restricted  to  specific  tactical  situations  or  sensor/weapon  systems 
and  required  rather  awkward  assumptions  about  the  nature  of  distribution  of  the 
various  attack  component  errors. 

The  effects  of  antisubmarine  weapons  as  a  function  of  distance  are  compara¬ 
tively  well  known.  Therefore,  the  major  difficulty  in  the  determination  of  kill 
probabilities  lies  in  the  calculation  of  the  probability  of  placing  a  weapon  within 
a  given  distance  from  the  target,  in  terms  of  the  expected  errors  in  the  various 
stages  of  attack. 

This  program  provides  a  distribution  of  these  distances  by  Monte  Carlo 
simulation  of  the  actual  antisubmarine  tracking  and  firing  tactics.  It  is  based 
on  a  model  (reference  (a))  which  is  both  flexible  enough  to  cover  most  tactical 
situations  and  sensor/weapon  systems  and  realistic  enough  to  yield  reliable  kill 
probability  values.  The  use  of  the  Monte  Carlo  technique  voids  the  need  for  un¬ 
reasonable  distortion  of  component  error  distributions. 


II.  GENERAL  DESCRIPTION 

Model 

The  tactical  situation  is  shown  in  figure  1.  The  attack  unit  is  at  the  origin 
of  the  coordinate  system.  The  true  course  of  the  submarine  is  parallel  to  the 
x-axis.  At  the  time  tj  *  0,  the  submarine  is  at  the  true  point  Pj  with  polar 

coordinates  (r^,  a)*.  Hie  submarine  proceeds  from  Pj  at  a  speed  Uj  and 

arrives  at  the  point  P2  at  time  t^.  It  then  executes  a  turn  of  radius  Rt  thru 

v  degrees  of  arc  at  a  velocity  U^.  When  the  turn  thru  y  is  completed  at  time 

tj  at  point  P^,  the  submarine  continues  on  a  tangential  course  at  speed  U^, 

reaching  point  P^  at  a  depth  %  at  time  t^ . 

The  time  t^  is  set  equal  to  the  actual  activation**  time  of  the  attack  weapon, 

thus  making  P^  a  variable  position  dependent  on  t^.  If  t^  is  less  than  t2, 

*'  AM  angles  are  measured  from  the  verttcal  or  the  range  line,  positive  direc¬ 
tion  to  the  right.  Deflections  arc  measured  from  the  range  line,  positive  to 
the  right. 

**  Activation  here  denotes  the  step  in  which  the  weapon  begins  to  exert  its  ef¬ 
fect  on  the  target  (c.g.,  beginning  of  search  for  a  homing  torpedo  or  detona¬ 
tion  for  a  depth  charge). 


P4  will  be  somewhere  between  Pj  and  Pj.  If  is  greater  than  t<2  but  less 
than  r^,  P^  will  be  somewhere  on  the  circular  arc  between  P,?  ami  P^. 

Finally,  if  t4  is  greater  than  t^,  P^  will  be  located  somewhere  along  the 
tangent  to  the  circular  arc  depending  on  the  time  elapsed  between  t^  and  t^. 

Due  to  errors  in  tracking,  the  attack  unit  presumes  the  submarine  to  be  at 
the  point  P'j  with  the  polar  coordinates  (r'j ,  <*')  at  time  t* ^  =  t  ^  =  0*.  The 

attack  unit  loses  contact  with  the  submarine  at  P'j  where  it  has  measured  the 

submarine's  speed  as  U'j  and  its  course  as  6'.  The  attacker  assumes  that 

the  submarine  continues  from  P'j  at  the  speed  U'j  on  the  course  6'  and 

calculates  P'j,  where  P'j  Is  the  assumed  location  of  the  submarine  at  t'j, 

the  desired  time  of  activation  of  the  weapon. 

The  aimpoint  of  the  weapon  P'g  is  specified  by  an  offset  angle  Ti  *  and 
an  offset  distance  h'  from  the  predicted  point  P'j.  1  ae  to  weapon  placement 
errors,  the  weapon  arrives  at  the  point  P^  at  depth  Z',  and  activates  at  time 
t'^.  The  time  t4  is  set  equal  to  t'4  so  that  the  distance  between  P^  and  P*4 
becomes  the  miss  distance  between  the  submarine  and  the  weapon. 

Program 

A  miss  distance  dj  is  calculated  for  each  Monte  Carlo  iteration.  At  the 
end  of  the  iteration  process,  two  graphs  are  plotted: 

Graph  1:  The  miss  distances  d(  are  plotted  in  sequential  ascending  order. 

Graph  2:  The  cumulative  frequency  distribution  of  the  miss  distances  is 
plotted  as  follows:  DM  is  an  input  to  the  program  and  is  used  to 

specify  the  largest  value  on  the  x-uxis.  If  Dj.  =  0  or  is  not  spec¬ 
ified,  the  largest  dj  or  dn  is  found. 

D..  is  then  chosen  as  the  smallest  number  of  the  form  1  x  10n, 

M 

2  x  10  ,  5  x  10,  -2  *  n  *  38,  which  is  greater  than  d,„.  .  Once 

UidA 

D..  has  been  chosen,  the  x-uxis  of  Graph  2  is  divided  into  N 

M  8 

equal  increments,  where  N,  is  also  an  input  to  the  program. 

b 

c  Primecf variables* are  used  to  describe  quantities  in  the  observed  system.  Un- 
primed  variables  denote  quantities  in  the  true  or  actual  system. 
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Each  increment  represents  a  range  of  miss  distances.  The  fre¬ 
quency  of  occurrence  for  each  range  of  miss  distance  is  calculated 
and  die  cumulative  frequencies  are  plotted. 

The  largest  and  smallest  values  of  d^  and  the  number  d^  >  are 

printed  out.  All  values  of  d^  >  are  ignored;  they  are  not  plotted 

on  Graph  1  and  they  are  not  used  in  determining  the  cumulative  fre¬ 
quency  distribution  of  Graph  2. 


III.  METHOD  OF  SOLUTION 


The  inputs  for  the  program  consist  of  the  following  parameters: 


Address 


Symbol 


Description 


1  ID  Identification  number  which  can  be  assigned  to  a  computer 

run  and  which  will  appear  on  the  graphs  and  the  printed 
output.  ID  is  an  integer  such  that  0  s  ID  s  32,  767 

2  Nj  Number  of  Monte  Carlo  iterations. 

3  Dw  The  value  of  miss  distance  that  will  be  used  as  the 

M  largest  value  on  die  x-axis  of  Graph  2  (feet)* 

4  Ng  Number  of  equal  increments  on  the  x-axis  of  the 

a  cumulative  frequency  plot  (one  increment  has  length  <* 

W*> 

5  N«  Number  of  "empty”  passes  thru  random  number 

generators.*** 

6  Tj  Distance  from  attack  unit  to  Pj,  die  submarine's  true 

position  when  the  last  contact  is  made  (feet) 

7  Atj  Observed  position  bias  error  in  range  (feet) 

8  or^  One  standard  deviation  of  position  error  in  range  (feet) 


If  not  entered,  TC  is  chosenas  die  smallest  number  of  the  form  1  x  10, 

2  x  10°,  5  x  10n  -2  <  n  s  38  which  is  greater  than  die  largest  calculated 
miss  distance. 


**  If  not  entered,  Ng  will  be  set  equal  to  500. 

***  Ng  passes  thru  the  random  number  generators  are  made  only  for  the  first 
data  set  processed  by  the  program. 
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Symbol 


Description 


*6 


a6 

U1 

u2 

u3 


>u 


1 


n* 


Radius  of  turn  that  submarine  executes  at  P2  (yards) 

True  bearing  from  attack  unit  to  submarine  at  point 
of  last  contact (degrees) 

Observed  position  bias  error  in  bearing  (degrees) 

One  standard  deviation  of  position  error  in  bearing 
(degrees) 

Number  of  degrees  of  arc  in  turn  of  radius  rt. 

If  positive,  o.,  is  one  standard  deviation  of  error  in 

the  turn  angle  v  (degrees).  >  will  be  replaced  by  a 
normal  distribution  about  .  If  negative,  v  will  be 
replaced  by  a  uniform  distribution  from  0  to  v. 

If  zero,  v  will  not  be  changed. 

Observed  course  bias  error,  where  true  course 
8  ■  90  (degrees) 

One  standard  deviation  of  course  error  (degrees) 

True  speed  of  the  submarine  from  Pj  to  P2  (knots) 

True  speed  of  the  submarine  from  P2  to  P3  (knots) 

True  speed  of  the  submarine  from  P^  to  P^  (knots) 

Observed  bias  error  of  submarine  speed  Uj  (knots) 

One  standard  deviation  of  error  in  submarine  speed 
Uj  (knots) 

True  time  at  which  submarine  starts  turn  (seconds) 

If  positive,  c  is  one  standard  deviation  of  error  in 

the  time  tj  (seconds),  tj  will  be  replaced  by  a  normal 

distribution  about  tj.  If  negative,  will  be  replaced 

by  a  uniform  distribution  from  0  to  to.  If  zero, 
tj  will  not  be  changed.  1 

Desired  time  of  activation  of  the  weapon  (seconds) 

Aimpoint  offset  angle  (degrees) 
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Address 

Symbol 

Description 

26 

h’ 

Aimpoint  offset  distance  (feet) 

27 

Weapon  placement  bias  error  in  range  (feet) 

28 

°v 

One  standard  deviation  of  weapon  placement  error  in 
range  (feet) 

29 

Aed» 

Weapon  placement  bias  error  in  lateral  displacement 
(feet)* 

30 

a 

ed* 

One  standard  deviation  of  weapon  placement  error  in 
lateral  displacement  (feet)* 

31 

Weapon  placement  bias  error  in  bearing  (degrees)* 

32 

o. 

One  standard  deviation  of  weapon  placement  error  in 
bearing  (degrees)* 

33 

Z 

True  depth  of  submarine  when  weapon  is  activated 
(feet) 

34 

AZ 

If  o  2  is  positive,  Az  is  considered  to  be  die  observed 

depth  bias  error  (feet).  If  oz  is  negative,  Az  is  con- 

sidered  to  be  the  maximum  operating  depth  of  the 
evading  submarine  (feet). 

35 

°Z 

If  positive,  o„  is  one  standard  deviation  of  error  in  the 

depth  Z  (feet).  The  observed  depth  of  the  submarine, 

Z',  will  be  calculated  as  a  normal  distribution  about 

Z  +  AZ  (feet).  If  negative,  Z*  will  be  calculated  as  a 
uniform  distribution  from  0  to  AZ,  the  maximum  operat¬ 
ing  depth  of  the  submarine  (feet).  If  zero,  Z'  will  be 
set  equal  to  Z  (feet). 

36 

V* 

Velocity  of  die  weapon  (feet/seconds). 

In  die  following  discussion,  is  a  Gaussian-distributed  (mean  *  0,  standard 

deviation  *  1)  random  number  which  is  always  less  than  4  standard  deviations. 

*u  is  a  uniformly-distributed  (mean  ■  0)  random  number. 


▼ 


if  ■  o„  *0, 
d  *d. 

If  icr  .  .  0, 


then  Option  2  (see  page 
then  Option  1  (see  page 


8)  will  be  used. 
8)  will  be  used. 


* 


* 

t 

ft 

A- 


t 


The  miss  distance  d  will  be  computed  once  the  x-y  coordinates  ol'  and 
P*4  are  determined. 

1.  It  is  first  necessary  to  calculate  the  x  and  y  coordinates  of  P^.  In 

order  to  do  this,  certain  variables  in  the  observed  system  must  be 
computed. 


Since  r{,  <5uj,  ?Uj,  Uj,  Uj,  and  u3  will  be  used  in  subsequent  calculations, 

they  must  be  scaled  to  be  consistent  with  the  units  of  the  rest  of  the  input 
parameters. 


1^*  1.6878065^,  i  *  1,  2,  3 
4Uj  -  1. 68780654  Uj 
0TI  *1.6878065o„ 

Uj  Uj 

*t  “  3rt 

Then  U'j  •  Uj  +  4U}  +RN 

o'  *  a  +  4a  +  Rm 

Nj  <* 

r>,.rI+ir,+RNj  ^ 


6' *90  +  A6  +  Rn  (6*90°  in  the  true  system) 

It  is  now  possible  to  calculate  die  x  and  y  coordinates  of  P’3>  x'3  and  y'3, 
from  die  geometry  of  figure  1. 


x'3  *  r'j  sinor*  +  U'jt'3  sin  6'  +  h'  sinT)' 
y'3  *  r'j  coso*  +  U'jt'g  cos  6'  +  h'  cos vf 

2.  The  next  step  is  the  determination  of  x'4  and  y'4,  the  x  and  y  co¬ 
ordinates  of  P*4«  This  calculation  depends  on  how  the  weapon  placement 
error  is  described: 

Option  1:  The  weapon  placement  error  is  described  by  range  and  deflection 
(sr,  and  -:d,  respectively) 

Option  2:  The  weapon  placement  error  is  described  by  range  and  bearing 
(ert  and  respectively) 
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Option  1 


V*‘V  +rn5 

ed's^ed*  +  RN  ce  (Cd*  *8  measured  positive  to  the  right  along 

^  ^  the  range  line  from  the  attack  unit  to  the 

submarine) 


rW(xV+<yV 

*'4  ■  *’3 + v  !!i +  'a-  Xs  * 


r4  */<xV  +  <*V 


Option  2 


V  -  AV  +  rn5  °€r, 
+  rn6  %. 


f3  •  ^3)  +  <7*3) 

r*4  *  r*3  +  V 
x*a 

/3*  »  arctan  — ^ 

^  3 

x*4*r*4  sin(/3*  +  €^,) 
y*4  *  r'4  cos  <0*  +  €g,) 


*  The  derivation  of  these  equations  is  somewhat  lengthy  and  is  therefore  not 
included. 


1  *  i 
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The  time  of  activation  of  the  weapon  is  calculated  taking  into  account  the 
delay  due  to  weapon  placement  error. 


3.  The  x  and  y  coordinates  of  P^,  and  y^,  must  now  be  determined: 

Since  tj  and  v  will  be  used  in  subsequent  calculations,  they  must  first  be 
transformed  in  accordance  with  the  convention  defined  on  page  5. 

U\  >0,  thenTj-tj+R^  o^ 

If  o  <  0,  then  T2  *  Ry  tj 
*2  1 

If  o^*0,  thenTj-tj 
IfoY>0,  thenr*Y  +  RNg 

If  o  <  0,  then  /*  *  Rtt  v 

'  u2 

If  ■  0,  then  r  ■  v 
Since  T2  and  r  have  now  been  datermined,  t^  is  given  by: 

«j*t2  +m  roj 

As  was  noted  previously,  the  location  of  P.  is  dependent  upon  t’ ;  thus,  there 
exist  three  possibilities: 

a)  f4  *  T2 

then  x4  *  rj  sin  a  +  Ujt'4 
y4  ■  Tj  cose 
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b)  T2  <  f4  *  t3 


then  let  9  =  (  *4  ~J2  )  JT 

v  l3  *2  y 

x4  » tj  sina  +  11^2+  Rt  sin  0 
y4  *  Tj  cosa  —  R{  (1  -  cos  0) 


c)f4>t3 

then  let  S  *  Ug  (t'4  —  tg) 

x4  *  rj  sina  +  UjTj  +  Rt  sin  y  +  S  cos  Y 
y4  *  cosa  -  Rt  (  1  -  cos  y )  —  S  sin  Y 


4.  The  final  step  is  the  calculation  of  the  miss  distance: 


If  oz  >  0,  then  Z'  -  2  +  AZ  +  Rj^  o2 


If  oz  <  0,  then  Z'  ■  Ru 
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AZ 


If  o2  ■  0,  then  Z*  *  Z 

Then  d  -  J(xA  -  x*4)2  +(y4  -  y'4)2  +(Z  -  Z')2 


4 


1 

\ 


I 


f 


4  See  page  6  to  review  the  definition  of  oz , 


1 
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IV.  USER’S  INSTRUCTIONS 


See  die  list  of  input  parameters  on  page  4.  The  user  should  submit  a  data 
submittal  form  (see  appendix  C).  These  forms  have  space  for  die  submitter's 
name,  the  date,  the  program  number  (13*64F),  the  security  classification  if  any, 
the  parameter  addresses  and  values,  and  an  estimate  of  the  computer  running 
time  (see  section  VIII). 

Input  flexibility  has  been  attained  by  allowing  the  user  to  vary  any  or  all  of 
the  parameters  in  a  computer  run.  In  any  one  set  of  data  the  parameters  remain 
fixed,  of  course,  but  there  is  no  programmed  limit  to  the  number  of  data  sets  a 
user  may  submit  in  a  run.  The  only  restriction  is  that  each  data  set  must  ter¬ 
minate  with  one  blank  card,  and  the  last  set  in  the  run  must  terminate  with  two 
blank  cards. 

A  further  advantage  enjoyed  by  the  user  is  that  for  each  data  set  afcer  die 
first,  he  need  submit  only  those  parameter  values  in  a  set  that  are  different  from 
those  in  the  previous  set.  *  This  is  accomplished  by  identifying  each  parameter 
by  a  '  nique  "address"  in  die  computer  memory  (see  appendix  E).  Initially  every 
address  is  cleared  to  zero  so  that  only  non  zero  input  parameters  need  be  entered. 


V.  SAMPLE  PROBLEM 

A  situation  was  fabricated  for  illustrative  purposes  whereby  most  of  the  bias 
errors  were  chosen  as  one  percent  of  the  true  value  and  standard  deviations  were 
taken  as  ten  percent  of  the  true  value.  The  output  is  shown  in  appendix  D;  and  the 
problem  submittal  form  in  appendix  C. 


VI.  KEYPUNCH  INSTRUCTIONS 

See  the  list  of  input  parameters  on  page  4  and  the  sample  problem  submittal 
form  in  appendix  C.  Input  card  formats  are  described  in  appendix  E. 


VII.  OPERATOR’S  INSTRUCTIONS 

Run  under  control  of  the  Bell  System  on  the  IBM  7090  computer.  Graphs  are 
made  by  die  Stromberg-Carlson  microfilm  recorder.  A  charactron  (CRT)  tape 
must  be  mounted  on  tape  unit  A8  in  the  low-density  mode.  When  die  program  is 
run  at  the  David  Taylor  Model  Basin  all  that  is  necessary  for  plotting  is  a  "BCRT" 
control^  card  (columns  8  thru  12)  and  the  "AMPLOS"  subroutine.  When  the  pro¬ 
gram  4crun  at  CEIR,  NAVIC  or  NAVCOSSACT,  die  "BCRT"  control  card  is  not 
to  be  used  and  the  "AMPLOS"  subroutine  is  to  be  replaced  by  the  "CXPLOT"  and 
"AMPLOF"  subroutines;  die  CRT  tape  on  A8  can  then  be  taken  to  DTMB  or 
Stromberg-Carlson  for  processing. 

*  The  exception  to  this  rule  is  N„,  the  number  of  empty  passes  through  the  ran¬ 
dom  number  generator  (see  pagb  4). 
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VIII.  TIMING 


The  table  presented  below  gives  a  conservative  estimate  of  the  computer 
running  time  for  one  data  set  as  a  function  of  N.,  the  number  of  Monte  Carlo 
iterations. 


(Number  of  Monte -Carlo 
Iterations) 


Computer  Running  Time 
(in  minutes) 


100  2 

1000  2 

2000  3 

3000  3 

4000  4 

5000  5 

6000  5 

7000  6 

8000  6 

9000  6 
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Reference:  (a)  Operations  Evaluation  Group,  "Design  of  Antisubmarine  Attack 
Models, "  in  preparation. 
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APPENDIX  B 

FORTRAN  STATEMENTS 


COMPUTER  PROGRAM  1  3-64P  ANT  I-SUbMAR  JNt:  ATTACK  AH/SAD 

DIMENSION  X 1(9000) *00(9000 > ,X2(  1 000  >  . VV (  1 000 >  ,00 (  36 )  .Gl  <  1 2  >  ,  62  (  S ) 
EQUIVALENCE  ( 00 ( 3) tOSUBM ) , ( 00(6) »RSUP1 ) » ( 00( 7  > ,DCLR1 ) , (00(H) ,S!Gwl 
X ) « (00(9 ) *RT  ARN) • (00( 1 0  >«  ALPHA ) * (0Q( 1 1 ) .OELALP) . (O0( I ? ) .SIGALP) , 
X(QO( 13) # GAMMA ) • ( 00( 14 ) , S 1 GGAM ) , (0Q( IS), OELDEL ) • ( 00 ( l(. ) .StGDEL ) . 
X(QO(  17)  ,U1KN0T  )  ,  (  0Q(  IS)  ,U2KN0T )  ,  (  00  (  19)  ,U3KN0T)  ,  {  CO(J’O)  ,  OELUA)  , 

X( 00(21 ) *S 1GUA )•( 00(22 ) *T IME2 )• (00(23) *SIGT2)*(Q0(?4)»T3P)* 

X(00(2S) *ETAP) • (00(26) • At  TCHP) • (00(27 )*0£L£RP) • (00(28) «  SIGERP) • 
X(Q0(29) *0£LEDP ) * (00( 30)* S1GE0P) • (Q0( 31 ) *OEL£BP) • ( 00 ( 3c  >  »  S l GEBP ) « 
X(00( 33) *ZEE ) • (00( 34) tOELZ ) • (QQ( 35) *S1GZ ) • (0O( 36) • VPRt  ME ) 

00  10  1*1*36 
10  00(I)«0. 

KEE*0 

NSUBS-900 

CALL  LGCHAR(8«4H  SAD) 

Pt*3* I41S9 

S  OK  121*740610306060 

5  41  ( 1 1 )*2 79  1  2147  3060 

a  G 1 ( 10 )*0 16060606060' 

8  «1(9)*6031 24601360 

8  GK  7 1*606060606060 

B  GK  6  >*443162626024 

B  Gl (9 >*3162632 14823 

8  Gl ( 4  >*256260314960 

8  Gl( 3 >*216223294924 

8  G 1 ( 2 ) *314527604681 

8  GK  I  >*242591606060 

0  G2( 19)*741006306060 

G2(14)*G1(U) 

8  G2( 13 ) *026060606060 

G2( 12 )*G1 (9 ) 

G2( 10)*G1 (7 ) 

0  G2( 9  >*23644464432 1 

8  G2( 8 )«633l 68296026 

6  G2( 7  >*512550642949 

8  G2( 6 >*2370602431 62 

a  62(91*635131226463 

8  62 ( 4  >*314645604626  . 

8  62 (3 >*6044 3 1626260 

a  G2( 2 >*243162632 148 

8  G2(! >-232562606060 

20  CALL  0ATA(Q0*1ND) 

1P( INO) 90*50 *30 
30  PRINT  40 
40  FORMAT (1H1) 

CALL  LGCHARI 8* 4H  SAD) 

CALL  ENDJ08 
90  IDENT*QQ( 1 ) 

NSU8I*QQ( 2) 

NSUBS*0Q(4> 

IF(KEE)  92*51*52 
51  NEMPTV*0Q(S) 

82  PRINT  60 

60  FORMAT  (1H1*  41H  CNA  PROGRAM  13-64P  ANTI-SUBMARINE  ATTACK  /// 

X8H  ADDRESS,  SX*  5H VALUE.  I  IX*  1 1 HOb SCRIPT I ON/) 

GEMPTV-NEMPTV 

PRINT  70* ( I «QQ( I >»I»1*4), GEMPTY* ( 1 ,001 1 ) , 1*6* 10) 

70  FORMAT ( IS*F 15*4*9X*  21HI0ENTIF ICATION  NUMBER/ 

XIS,F15*4*5X« 29HN0*  OF  MONTE-CARLO  ITERATIONS/ 
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XI5, FIS. A.SX  .35HLARGEST  VALUE  FOR  X-AXIS  ON  GRAPH  ?./ 

XI5. FIS. A.SX .37HNUMUER  OF  EOUAL  INCREMENTS  ON  GRAPH  it 
XSH  5,FI*i.A.5X,33HMUMBF.R  OF  EMPTY  RANDOM  NO.  PASSES/ 
XI5.FIS.4. 5X  ,  7HH  SUB  1/ 

XI5.FIS.A.SX.I3HOELTA  R  SUB  1/ 

XI5.F IS.A.SX. l JHSIGMA  R  SUH  1/ 

XIS.F1S.4. 5X . 7HR  SUB  Tt 
XI5.F1S.A.SX »SH ALPHA/ ) 

PRINT  80.11 .0Q( t). 1*11,20) 

80  FORMAT! IS.F IS.A.SX. 1 IhOELTA  ALPHA/ 

XIS.FIS.A.SX, 1SHS1GMA  SUH  ALPHA/ 

XIS.FIS.A.SX, SHGAMMA/ 

X15.F  IS.A.SX  .  ISHSIGMA  SUli  GAMMA/ 

XI3.F IS.A.SX. 1 IHDCL  T  A  DELTA/ 

XIS.F IS. A. SX . ISHSIGMA  SUH  DELTA/  ' 

xts.r is.a.sx ,7hu  sub  i/ 

XIS.FIS.A.SX.7HU  SUB  2/ 

XIS.FIS.A.SX. 7HU  SUB  1/ 

XIS.FIS.A.SX.t  JHUELTA  U  SUB  I/) 

PRINT  90. I  I .001 | ) , 1*21,30) 

90  FOHMAT  IIS. FIS. A.  SX. I  7HSI  GM,\  SUH  U  v»U0  1/ 

XIS.FIS.A.SX. 7HT  SUB  it 

XIS.F IS. A.SX . | 7HSIGMA  SUB  T  SUB  2/ 

X  IS »F IS.A.SX . I JHT  PRIME  SUB  3/ 

XIS.FIS.A.SX, 9HF:TA  PRIME/ 

XIS.FIS.A.SX, 7HH  PRIME/ 

XIS.FIS.A.SX .23HPELTA  EPSILON  SUH  R  PRIME/ 

XIS.FIS.A.SX .29HS IGMA  SUB  EPSILON  SUH  R  PRIME/ 

XIS.F IS.A.SX .2SHDELTA  EPSILON  SUH  0  PRIME/ 

XIS.F IS. A »SX .29HSIGMA  SUB  EPSILON  SUH  O  PRIME/) 

PRINT  100,1 1.001 11.1*31.36) 

100  FORMAT!  IS.F- lb.  A. 5X.2BH0ELTA  EPSILON  SUH  BETA  PR  IMF/ 
XIS.FIS.A.SX.32HSIGMA  SUB  EPSILON  SUP  BCTA  PRIME/ 

xts.r is.a.sx. ih it 
XIS.FIS.A.SX. 7HDCL T A  It 
XIS. FIS. A.SX. IIHSIGMA  SUH  It 
XIS.FIS.A.SX ,7HV  PRIME////) 

IFINEMPTY)  liU.IA0.ll9 

110  00  120  l*|.  NEMPTY 
CALL  RANUMHIOUMMY ) 

120  CALL  GRNUMA IOUMMY) 

130  NEMPTVsO 

IAO  UIFPS* I ,687(J0f>S*U  1KN0T 
U2FPS* I <68/ A06S*U2KN0T 
U3FPS» 1 .68/80/ S«U3KN0T 
DELUIal .687ROOb«OELUA 
SI  GUI* 1 .6B/806S*SIGUA 
RTURN«3.*RT  AHN 
1)0  330  1st,  NSUBI 
CALL  GRNUMA ( RN I ) 

CALL  GRNUMA (RN2) 

CALL  GRNUMA I RN3) 

CALL  GRNUMA (RNA) 

UIP*UIFPS+DELU1FRNI*SIGUI 
ALPHAPs  ALPH AFDEL ALPFRN2*S igalp 
RIPsRSUBt FOELRI 4RN3*S IGR l 
OELP*90 •♦DtLOEL  4RNA4S I GDEL 

X3P*R IP4SINBF I ALPHAP )*U|P  »T3P»S I NDF  <  DFLP ) ♦ AI T  CHP»S INOF (ETAP 
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V3P*RIP9C0S0F! ALPHAPH'UIPtT3P*C0SDF<0ELP)4AITCHP*C0SDF<ETAP) 

CALL  GRNUN4  4  RNS ) 

CALL  GRNUN4 ( RN6) 

ERP*0ELCRP4RNS«S I GERP 
R3P* SORTP 1 X 3P4# 2 ♦ V 3P« # 2 ) 

IP  (OELEBPI  »?0. ISO. 170 
ISO  IP  ISIGEBP)  170. 160* iro 
CALCULATIONS  POR  OPTION  I 
ISO  E0P*0ELCDP*RN6*SIGE0P 

X4P»X3P4ERP*X 3P/R3P*E0P*Y3P/R3P 
Y4P«Y3P4ERP*Y3P/R3P"*EDP*X3P/R3P 
R4P*S0RTP(X4P«*24Y4P**2 ) 

SO  TO  ISO 

CALCULATIONS  POR  OPTION  2 
170  EBP*0ELE8P4  RN6*S 1 6EBP 
R4P«R3P4ERP 
ANGLE* ATANOP ( V  JP/X  3P  > 

IP  IXSP)  172.174.174 
172  ANGLE* ANGLE  +  180 • 

174  BCTAP*90.-ANGLE 

X4P*R4P*SIN0P( BETAPtEBP) 

V4P*R4P*C0SDF! BETAPAEflP) 

180  T4P*T3P*!R4P-M3P)/VPRIME 
TTf«TIME2 

IP  (SIGT2)  t«0. 210. 200 
190  CALL  RANUMO  (RUI) 

TT2*RU  I  *T IMC2 
GO  TO  210 

200  CALL  0RNUN4 ( RN7 ) 

TT2*T INE24RN7*SIGT2 
210  GGANMA*GANMA 

IP  (SIGGAMI  220.240.230 
220  CALL  RANUNB  IMU2) 

GGAPNA«RU2*GANMA 
GO  TO  240 

230  CALL  GRNUN4IRN8) 

GGANNAaGANM A+RNH*S IGGAM 
240  IP  (T4P-TT2)  2S0.2S0.260 
280  XPOURaRSUBt  OSINOPI ALPHA) ♦U1FPS»T4P 
VP0UR«RSUB1 4C0S0P! ALPHA) 

60  TO  290 

260  TIME3«TT24(P|/inO.)*r.GAMMA*RTURN/U2FPS 
1PIT4P-TINE3)  270.270.2fl0 
270  THETA*! (T4P-TT2)/! TIN63-TT2) )*GGAMHA 

XP0UR*RSU81 «SI NOP (ALPHA ) ♦UIFPS»TT24RTURN*S1NDF( THPTA) 

VPOURaRSUBt  *C0Sl>P  I  ALPHA ) -RTURN*  (  I  .-COSOP!  THETA  )  ) 

GO  TO  290 

280  ESS*U3PPS*!  T4P-TINE3) 

XPOURaRSUB I • S I  NOP  I ALPHA ) tU IPPS»T T2*R TURN*  S I NOP ! GG ANNA )4ESS*C0S0P 
XI6GANNA) 

VPOURaRSum *00307 ( ALPHA ) -RTURN*! t .-COSOP ( GGAMMA ) )-ESS*SINOP 
X1GGANNA) 

290  ZPRINE*2EE 

IP  (SIG2)  300.320.310 
300  CALL  RANUN8IRU3) 

2PH!ME*RU3«OEL< 

GO  TO  320 

310  CALL  GRNUN4IRN9) 
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**** 


ZPHIME*ZEE40FLZ4«N9»SIG2 

320  00< I )*SORTF  ( < XF0UP-X4P >**24( YFOUR-V4P > **24 ( ZKE-ZPRIMC )**2> 

330  XI ( I )* I 

CALL  FORAST (OD. NSUBI »  l  ) 

OMAX*CO ( NSUB I ) 

SCALE*0SU8M 

IF  (OSUBM)  340 . 340*360 
340  00  370  K*l»  41 

SCALE** 00 l* 10*  **K 
IF  (OMAX-SCALE)  380.380.350 
350  SCALE*2.*SCALE 

IF  { OMAX-SC ALE )  380.380.360 
360  SCALE*? »S*SC ALE 

IF  (OMAX-SC ALE)  380*380.370 
370  CONTINUE 
380  J 33*0 

GGG«NSUBS 
OELX*SCALE/GGG 
00  390  1*1. NSUBS 
390  YVII)*0. 

00  420  1*1. NSUBI 
INl)EX*OD(  1 1  /0ELX40.99999999 
IF  (INOEX-NSURS)  410.410.400 
400  333*33341 
GO  TO  420 

410  YV ( INDEX )*V V( I NPEX 1 4 1 . 

420  CONTINUE 

GNORM*NSUB 1-333 
X2(l)*0bLX/2. 

00  430  3*2. NSUBS 
X2 (3)*X2(3— 1) 40ELX 
YV (3)*YV(3) 4YY  4  3-1) 

430  YV!  3- 11*100*  *VY( 3— 1 1 /GNORM 

YY ( NSUBS ) *1 00* • YY ( NSUBS ) /GNORM 
GI(8)«BIN0F| IUbNT.6) 

G2( 1 1 )*G1 (8 ) 

XUPPER*Xl (NSUBI 1 
I NORM*GNORM 

CALL  FNPLOT ( G1 ( 121  * 12H( 8H  NUMBER) . ?9H( 23H  MISS  DISTANCE  IN  FEET 
X>. 1..XUPPER.0*. SCALE. 20. 20.20.1 .6MIF5.0) .6HIF6.0) ) 

CALL  CURVE! XI ( I  NORM) *0D( (NORM) , IN0RM.6H  ) 

CALL  FNPLOT <G2( 15). 28H(?3H  MISS  DISTANCE  IN  FtET)»l3H(9H  PERCENT 
X) .0*. SCALE. 0*. 100. .20* 2* 20. 1.6H(F6.0)*6H(F4.0>) 

CALL  CURVE! X2INSUBS) • YY(NSUBS) .NSUUS.6H  ) 

PRINT  440.0MAX. 00! 0.333. SCALE 

440  FORMAT (29H  THE  MAXIMUM  MISS  OISTANCE  •  F9.2.5H  F6F.T/29H  THE  MINIMU 
XM  MISS  OISTANCE  *  F9.2.SH  FEET//12M  THERE  MERE  I4.30H  MISS  OISTAN 
XCES  GREATER  THAN  F9.2.SH  FRET) 

KEE*  l 
60  TO  ?0 
END 
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APPENDIX  C 

SAMPLE  PROBLEM  SUBMITTAL  FORM 


OBO  COMPUTER  DATA  8UMITTAL  FORM 

_ Data:  !?<>*/ 

Program  No.  I~?>  -  Hi'P  Bat.  Tima  3  Claaalflcatioa  (I  Aft  I 

jpjdti  I— ciulUoih  ■■  "" 


Addraoa 

Valua 

Addraaa 

Valua 

Addreaa 

Valua 

Addraaa 

Valua 

1 

(S*o 

wm 

3 

A 

1*00 

ra 

/t* 

0 

wm 

/« 

IB 

4* 

J*6 

wm 

3“ 

IBP 

ira 

r _ 

rat 

... 

WM 

4 

ph mi 

wm 

S* 

1 

/» _ 

wm 

/• 

i r 

PS9H! 

wm 

J* 

1 

_ 

wm 

4 

H 

m 

j* 

wm 

6 

mm 

■X 

wm 

3 

wm 

X  ■ 

3X 

m 

_ 

4* 

33 

warn 

wm 

4 

\rm 

.3 

wm 

f 

\wm 

31 

wm 

9 

i  rm 

4* 

warn 

J9 

— — 

Ib 

wm 

J* 

^  — 

wm 

3 * 

m 

•3 

NOTES: 

1.  A  valua  of  aaro  muat  ba  oatarad  u  0.  aotlaftblaak. 

*.  Dacinal  pta.  may  ba  oolttad  U  uadaratood  to  follow  tha  rtpmoat  dpt. 

9.  Tha  yah*  3  X  10'5  may  ba  aatarad  aa  .00003  or  3-9,  aotaa  9  X  10  . 

4.  Tha  factor  portion  of  a  valua  may  aot  eoatala  a>or«  thaa  I  dlgiu. 

3.  Tha  aapoaaat  portloa  of  a  valua  moat  Ua  arttkia  lha  raagt  an. 

6.  Bapooaata  may  ba  oatttad  H  aaro.  V  aot.  (bay  muat  ba  aigaad. 

7.  Haak  carda  abould  ba  ladicatad  by:  —  b . ■  C-l 

_ _  (REVERSE  BLANK) 
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APPENDIX  D 

SAMPLE  PROBLEM  OUTPUT 

CNA  PROGRAM  I3-64P  ANTI-SUBMARINE  ATTACK 


ADDRESS 

VALUE 

DESCRIPTION 

* 

1 

1500.0000 

IDENTIFICATION  NUMBER 

2 

2000.0000 

NO.  OF  MONTE-CARLO 

ITERATIONS 

3 

0. 

LARGEST  VALUE  FOR 

X-AXIS  ON  GRAPH  2 

4 

100.0000 

NUMBER  OF  EQUAL 

INCREMENTS  ON  GRAPH  2 

5 

5.0000 

NUMBER  OF  EMPTY 

RANOOM  NO.  PASSES 

6 

5000.0000 

R  SUB  1 

7 

50.0000 

OELTA  R  SUB  1 

8 

500.0000 

SIGMA  R  SUB  1 

* 

9 

100.0000 

R  SUB  T 

10 

20.0000 

ALPHA 

11 

0.2000 

OELTA  ALPHA 

12 

2.0000 

SIGMA  SUB  ALPHA 

13 

60.0000 

GAMMA 

14 

6.0000 

SIGMA  SUB  GAMMA 

* 

IS 

0.9000 

DELTA  OELTA 

16 

9.0000 

SIGMA  SUB  DELTA 

17 

30.0000 

U  SUB  1 

16 

20.0000 

U  SUB  2 

19 

30.0000 

U  SUB  3 

20 

0.3000 

OELTA  U  SUB  I 

21 

3.0000 

SIGMA  SUB  U  SUB 

1 

22 

100.0000 

.  T  SUB  2 

23 

10.0000 

SIGMA  SUB  T  SUB 

2 

24 

300.0000 

T  PRIME  SUB  3 

2S 

45.0000 

ETA  PRIME 

26 

50.0000 

H  PRIME 

27 

10.0000 

DELTA  EPSILON  SUB 

R  PRIME 

26 

50.0000 

SIGMA  SUB  EPSILON 

SUB  R  PRIME 

29 

0. 

OELTA  EPSILON  SUB 

D  PRIME 

30 

0. 

SIGMA  SUB  EPSILON 

SUB  0  PRIME 

31 

3.0000 

OELTA  EPSILON  SUB 

BETA  PRIME 

32 

2.0000 

SIGMA  SUB  EPSILON 

SUB  BETA  PRIME 

33 

50.0000 

Z 

34 

0.5000 

DELTA  Z 

35 

5.0000 

SIGMA  SUB  Z 

36 

600.0000 

V  PRIME 

THE  MAXIMUM  MISS  DISTANCE  *  17726.80  MEET 
THE  MINIMUM  MISS  DISTANCE  ■  1012.03  FEET 

THERE  WERE  0  MISS  DISTANCES  GREATER  THAN  20000.00  FEET 
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APPENDIX  E 


DATA  SUBROUTINE 


1 .  Introduction: 

Many  computer  programs  require  the  flexibility  of  varying  any  or  all  of  the— 
parameters  in  a  computer  run.  Although  FORTRAN  is  fairly  flexible  in  its 
arithmetic  and  control  statements,  its  input-output  statements  are  quite  rigid. 

In  order  to  read  cards  for  instance,  considerable  effort  must  be  expended  by  the 
FORTRAN  programmer  in  writing  his  input  statements.  This  subroutine  eliminates 
some  of  that  tedium.  The  concept  of  a  "data  set"  is  used.  A  data  set  consists  of  a 
sequence  of  punched  cards  terminated  by  one  blank  card.  A  parameter  deck  for  a 
computer  run  may  consist  of  several  data  sets.  Such  a  parameter  deck  is  terminated 
by  two  blank  cards. 

2.  Parameter  Addresses: 

The  primary  advantage  of  this  subroutine  over  FORTRAN  input  statements  re¬ 
sults  from  the  use  of  "parameter  addresses."  An  address  is  a  relative  location  in 
the  computer  memory.  It  is  the  subscript  of  an  array  or  matrix.  For  example,  in 
an  array  called  X  ,  the  parameter  value  X^  would  be  located  at  address  53 .  By 

using  the  parameter  addresses,  a  user  of  the  program  need  submit  only  those  param¬ 
eter  values  in  a  data  set  that  are  different  from  those  in  the  previous  set. 

Three  types  of  addresses  are  permitted  by  this  subroutine. 

(1)  A  numeric  address  consisting  of  one  to  five  characters,  each  of  which 
is  a  digit  0*9.  Such  an  address  (n)  refers  to  the  nth  element  in  a 
specified  array. 

(2)  An  alpha  address  consisting  of  one  to  six  characters,  the  first  of  which 
must  be  alphabetic  (A-Z).  The  remaining  may  be  alphabetic  or  numeric 
(A-Z  or  0-9).  Such  an  address  refers  to  the  n“*  element  in  a  specified 
array  (Unf  26),  where  the  first  character  of  the  address  corresponds 
to  n  as  the  26  letters  of  the  alphabet  correspond  to  the  integers  1-26. 

(3)  A  matrix  address  consisting  of  two  or  more  numeric  fields  separated  by 
commas.  For  example,  the  address  53,  47  refers  to  the  element  in  the 
53fd  row  and  the  47t*»  column  of  a  two-dimensional  matrix.  There  is  no 
limit  to  the  number  of  dimensions  in  a  matrix  address. 

3.  Input  Card  Format: 

A  standard  submittal  form  (see  attachment)  has  been  designed  for  the  analyst. 
This  form  provides  for  entering  parameter  values  with  their  associated  addresses. 
The  user  indicates  blank  cards  to  separate  data  sets.  The  keypunch  opera  or  has 
the  option  of  punching  one  address  and  value  per  card,  or,  if  the  addresses  are 
sequential,  of  punching  one  address  and  several  values  on  a  card. 


Only  columns  1-72  of  a  card  are  used.  Etch  column  must  contain  one  of  the 
following:  a  digit  (0. 9),  a  "+"  or  "  -  "  sign  or  a  dash,  a  letter  (A-Z),  a  period, 
a  comma,  or  a  blank.  Each  punched  card  must  contain  one  parameter  address. 

The  address  may  start  in  column  1,  or,  if  desired,  may  start  in  a  later  column, 
provided  all  columns  before  it  are  blank.  The  address  is  terminated  by  at  least 
one  blank  column.  Only  one  address  is  permitted  on  the  card.  Succeeding  col¬ 
umns  contain  one  or  more  parameter  values,  each  separated  by  one  or  more 
blank  columns.  A  value  may  be  signed  or  unsigned.  The  length  of  the  value  field 
is  variable.  No  blanks  are  permitted  within  a  value  field.  A  value  may  be  punched 
with  or  without  an  exponent.  An  exponent  is  recognized  by  the  presence  of  a  plus 
or  minus  sign  (or  dash)  between  the  fractional  part  and  exponent  part  of  the  value. 
Decimal  points  (periods)  may  be  punched  in  either  the  fractional  or  exponent  parts 
of  a  value.  If  more  than  one  value  is  punched  on  a  card,  those  after  die  first  will 
be  entered  at  sequential  addresses  relative  to  the  address  of  the  first  value. 

4.  Usage: 

A  data  set  is  read  by  the  use  of  the  statement: 

CALL  DATA  (X,  I) 

in  a  FORTRAN  program  for  the  IBM  7090.  The  argument  X  is  the  name  of  an 
array  ir.  the  program.  The  argument  I  is  an  indicator  set  by  the  subroutine. 

This  indicator  may  be  tested  by  the  main  program  upon  return  from  the  sub¬ 
routine.  It  will  have  a  value  of  0  or  1  or  2 . 

0:  The  subroutine  has  read  a  data  set.  The  main  program  will  normally 
proceed  to  operate  on  this  data. 

1:  The  subroutine  has  read  the  second  blank  card  which  terminates  the 
parameter  deck.  The  main  program  will  normally  terminate  at  this 
point. 

2:  The  subroutine  has  read  a  "bad"  data  card.  The  main  program  may 
terminate  the  run,  or  ignore  the  card  and  return  to  the  subroutine  to 
read  the  rest  of  the  data  set. 

If  the  cards  to  be  read  contain  matrix  addresses,  additional  arguments  must  be 
included  in  the  FORTRAN  calling  statement: 

CALL  DATA  (X,  Dy  Dg,  Dg . Dr,  I) 

where  Di  is  the  ith  dimension  of  the  matrix  X. 

5.  Method: 

See  the  attached  flow  chart.  DATA  reads  parameter  values  and  loading  addres¬ 
ses  from  cards.  If  sense  switch  S  is  up,  it  will  read  the  values  and  addresses  from 
tape  (unit  A2).  It  converts  die  values  to  floating  point  numbers,  and  stores  them  as 
elements  of  an  array  specified  in  die  calling  statement.  The' elements  are  specified 
by  the  addresses.  If  a  card  (or  tape  record)  is  read  which  contains  non- permitted 
characters  (see  input  card  format  above),  DATA  prints  the  statement  "bad  data  card, 
followed  by  an  image  of  the  card  itself. 


6.  Coding  Information: 

See  the  symbolic  listing  included  in  this  appendix.  DATA  is  written  in  the  7090 
FAP  language.  It  must  be  used  in  conjunction  with  the  BELL  system.  It  requires 
401  words  storage  space. 


E-3 


• 


3naarr,T.H3.. 


TRA 

BAD 

LXA 

FIELD. 2 

TRA 

Fl+1.2 

TRA 

FT 

FIELD"?  (EXPONENT  FIELD) 

TRA 

F6 

FIELD-6  (VALUE  FIELD) 

TRA 

F5 

FIELD"*  (BLANKS  AFTER  ADDRESS) 

TRA 

F4 

FIELD"4  (ARRAY  ADDRESS) 

TRA 

F3 

F2£LD>3  (ALPHA  ADDRESS) 

TRA 

F2 

FIELD-2  (NUMERIC  ADDRESS) 

FI 

TXH 

8.4. 41 

F!ELD«1  (BLANKS  BEFORE  ADDRESS) 

F1A 

STZ 

•• 

RETURN  INDICATOR  >  0 

AXT 

2.2 

FIELD  •  2 

SXA 

FIELD. 2 

TXH 

H.4.J1 

•  NUMERIC  CHARACTER 

TXH 

1.4.28 

SIGN  OR  DASH 

TXL 

BAD .4 .2 

AXT 

3.2 

ALPHA  CHARACTER.  FIELD  >  3 

SXA 

FIELD. 2 

TXI 

•♦1 .4.-2 

ADDRESS  >  NTH  ALPHA 

SXA 

ADDRES.4 

TRA 

B 

BLANK  CHARACTER 

F2 

TXH 

J.4.41 

TXH 

H.4.31 

NUMERIC  CHARACTER 

TXH 

K.4.28 

SIGN  OR  DASH 

TXH 

BAD .4. 2 

TXH 

L.4.1 

PERIOD 

AXT 

4*2 

COMMA.FIELD  >  4 

SXA 

FIELD.2 

AXT 

2*2 

SXA 

DIMENS*2 

DIMENSION  ■  2 

cla 

ADORES 

TEST  ADDRESS 

TIE 

BAD  . 

TMI 

BAD 

F2A 

STZ 

ADDER 

AO0ER"O 

TRA 

B 

FS 

TXH 

J.4,41 

BLANK  CHARACTER 

TXH 

B.4.31 

NUMERIC  CHARAC’ER 

TXH 

K.4.28 

SIGN  OR  DASH 

TXH 

8.4.2 

ALPHA  CHARACTER 

TXH 

L.4.1 

PERIOD 

TRA 

BAD 

F4 

TXH 

M.4.41 

BLANK  CHARACTER 

TXH 

N.4.31 

NUMERIC  CHARACTER 

TXH 

P.4,28 

SIGN  OR  DASH 

TXH 

BAD .4, 2 

TXH 

0.4,1 

PERIOD 

TRA 

T 

COMMA 

F5 

TXH 

B.4.41 

BLANK  CHARACTER 

AXT 

6,2 

FIELD  •  6 

SXA 

FIELD.2 

TXH 

U.4,31 

NUMERIC  CHARACTER 

TXH 

K.4.28 

SIGN  OR  DASH 

TXH 

BAD, 4. 2 

TXH 

0.4,1 

PERIOD 
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TRA 

TXH 

BAO 

X»4*41 

BLANK  CHARACTER 

TXH 

U*4*31 

NUMERIC  CHARACTER 

TXH 

Y.4.2B 

SIGN  OR  DASH 

TXH 

TXH 

BAD** *2 
Gt**l 

PERIOD 

TRA 

TXH 

BAO 

2*4*41 

BLANK  CHARACTER 

TXH 

BB*4*31 

NUMERIC  CHARACTER 

TXH 

EE»4*2I 

SIGN  OR  DASH 

TXH 

TXL 

AXT 

BAD *4 *2 

BAO *4*1 

1*2 

PERIOD*  POINT  •  1 

SXA 

TRA 

LOO 

POINT  *2 

B 

ADORES 

ADDRESS  •  10  X  ADORESS  4  N 

XCA 

ACL  CHARAC 
STO  ADORES 

TRA  6 


TXH  6*4*30 

♦  SIGN 

CLA  ADORES 

SSM 

STO  ADORES 

TRA  6 

SET  SIGN  OF  ADDRESS 

AXT  3*2 

SXA  PIEL0*2 

TRA  B 

FIELD  ■  3 

TXH  LI *4* SO 

4  SIGN 

CLA  VALUE 

SSM 

STO  VALUE 

TRA  LI 

SET  SIGN  OF  VALUE 

AXT  1*2 

SXA  POINT *2 

POINT  •  l 

AXT  6.2 

SXA  FIELD.2 

TRA  B 

FIELD  •  6 

AXT  3*2 

SXA  FIELD.2 

TRA  S 

FIELD  •  3 

LDQ  ADDER 

MPV  HIO 

STO  ADDER 

TSX  T1.4 

MPY  CHARAC 

XCA 

ADO  ADDER 

STO  ADDER 

TRA  B 

ADDER  ■  10  X  ADDER  ♦  N  X  PROD 

TXH  R»4»30 

4  SIGN 

CLA  VALUE 

SSM 

SET  SIGN  OF  VALUE 

STO 

VALUE 

TRA 

R 

0 

AXT 

1.2 

SXA 

POINT .2 

R 

AXT 

6*2 

SXA 

FIELD. 2 

S 

LXA 

EXIT. 2 

TXI 

•♦1.2. -3 

PXA 

0.2 

SUB 

D I MENS 

TNZ 

BAD 

T 

TSX 

T1.4 

CLA 

ADDER 

SUB 

PROD 

STO 

ADDER 

TIE 

BAD 

TMI 

BAD 

ADD 

ADDRES 

STO 

ADORES 

CLA 

DIMENS 

ADD  HI 

STO  DIMENS 

TRA 

F2A 

T1 

SXA 

T4t4 

CLA 

HI 

STO 

PROD 

STA 

T3 

LXA  DIMENS*2 

TXI 

•♦1.2.-1 

LXA 

XA.4 

T2 

CAL 

T3 

ADD  HI 

STA 

T3 

T3 

CLA 

**.4 

STA 

•♦1 

LDQ 

ft* 

ROL 

IB 

MPY 

PROD 

STO 

PROD 

TIX 

T2.2.1 

TA 

AXT 

•*.4 

TRA 

1*4 

U 

CLA 

POINT 

tne  v 

LDQ  VALUE 

FMP 

DEC10 

SSP 

FAD  NUMB 

LDQ  VALUE 

LLS 

0 

STO  VALUE 
TRA  B 

V 

LXA 

POINT .4 

CLA 

NUMB 

POINT  «  1 
FIELD  «  6 
CHECK  DIMENSION 


ADDER* ADDER-PROD 

CHECK  ADDER 


PROD  ■  PRODUCT  OF  DIMENSIONS 


TEST  POINT 

VALUE  •  10  X  VALUE  ♦  N 


VALUE  ■  VALUE  ♦  N/(10*«POINT ! 
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FDP  DEC 10 
XCA 

TIX  *-2,4,1 
LOO  VALUE 
LLS  0 
FAD  VALUE 
STO  VALUE 
TRA  DO 

TXH  B.4.30  4  SIGN 

CLA  VALUE  SET  SIGN  OF  VALUE 

SSM 

STO  VALUE 
TRA  B 

CLA  XLOC  XI ADDRESS!  •  VALUE 

SUB  ADORES 

STA  *42 

CLA  VALUE 

STO  « 

TRA  AA 
AXT  7*2 

SXA  FIELD«2  FIELD  •  2 

TXH  00,4,30  4  SIGN 

CLA  EXR  SET  SIGN  OF  EXP 

SSM 

STO  EXP 
TRA  GG 

CLA  XLOC  XI ADDRESS)  •  VALUE  X  10**EXP 

SUB  ADORES 

STA  21 

CLA  DEC10 

LDO  EXP 

CALL  EXP (3 

XCA 

FMP  VALUE 
STO  ** 

ST2  VALUE  VALUE  •  0 

TRA  FF 

CLA  POINT  TEST  POINT 

TN2  CC 

LOO  EXP  EXP  •  10  X  EXP  4  N 

FMP  DEC10 

SSP 

FAD  NUMB 
LDO  EXP 
LLS  0  . 

STO  EXP . 

TRA  B 

LXA  POINT t*  EXP  »  EXP  ♦  N/l 10**PO!NT) 

CLA  NUMB 
FDP  DEC10 
XCA 

TIX  •«2»A»1 
LDO  EXP 
LLS  0 


E-ll 


FAD  EXP 
STO  EXP 

00  CIA  POINT 

AOO  HI 
STO  POINT 
TRA  B 

EE  CLA  XLOC 

SUB  ADORES 
STA  EE1 
CLA  0EC10 
LOO  EXP 
CALL  EXPO 
XCA 

FMP  VALUE 

EE1  STO  ** 

PXO  0.0 

TXH  *♦2*4.30 

SSH 

STO  VALUE 

FF  STZ  EXP 

AXT  5*2 
SXA  FIELD. 2 
CAL  ADORES 
ADO  HI 
SLW  ADORES 

GG  STZ  POINT 

TRA  6 

HH  LXA  FIELD. 1 

TXL  JJ.1.1 
TXL  BAD. 1.4 
TXL  A.1.5 
TXH  11.1,6 
CLA  XLOC 
SUB  ADORES 
STA  *+2 
CLA  VALUE 
STO  ** 

TRA  A 

II  TXH  BAO.1,7 
CLA  XLOC 
SUB  ADORES 
STA  III 
CLA  DEC 10 
LOQ  EXP 
CALL  EXPO 


III 

XCA 

FMP 

STO 

VALUE 

•• 

BAD 

TRA 

TSX 

A 

HPRINT.4 

112 

P2E 

AXT 

SXD 

PRINT, 0.15 

2.1 

**.i 

XI 

AXT 

**.i 

POINT  ■  POINT  ♦ 


1 


X( ADDRESS)  •  VALUE  X  10**EXP 


VALUE  ■  0 
♦  SIGN 

SET  SIGN  OF  VALUE 

EXP  ■  0 
FIELD  >  5 

ADDRESS  >  ADDRESS  ♦  1 


POINT  b  0 


FIELDbI,  EXIT 

FIELD-5.  READ  ANOTHER  CARD 
FIELDb6.  X( ADDRESS)  s  VALUE 


FIELD-7. 

XI ADDRESS)  -  VALUE  X  10**EXP 


E-12 
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X2  AXT  ft*»2 

X*  AXT  *ft*4 

EXIT  TRA  ft*»4 

MASK  OCT  777777TOOOOO 

PRINT  BCD  3  BAD  DATA  CARD.*. 


CARD  BSS 
ADORES  HTR 
VALUE  HTR 
EXP  HTR 
POINT  HTR 
PIELD  HTR 
COLUMN  HTR 
TABLE  OCT 

12 

•• 

H 

ft* 

*• 

*• 

ft* 

60 

BLANK 

OCT 

0 

0 

OCT 

1 

1 

OCT 

2 

2 

OCT 

3 

3 

OCT 

A 

4 

OCT 

5 

3 

OCT 

6 

6 

OCT 

7 

7 

OCT 

10 

8 

OCT 

11 

9 

OCT 

20 

♦  SION 

OCT 

40 

-  SION 

OCT 

14 

DASH 

OCT 

71 

Z 

OCT 

70 

Y 

OCT 

67 

X 

OCT 

66 

W 

OCT 

63 

V 

OCT 

64 

U 

OCT 

63 

T 

OCT 

62 

s 

OCT 

31 

R 

OCT 

30 

0 

OCT 

47 

P 

OCT 

46 

0 

OCT 

43 

N 

OCT 

44 

M 

OCT 

43 

L 

OCT 

42 

K 

OCT 

41 

J 

OCT 

31 

I 

OCT 

30 

H 

OCT 

27 

0 

OCT 

26 

F 

OCT 

23 

E 

OCT 

24 

D 

OCT 

23 

C 

OCT 

22 

B 

OCT 

21 

A 

OCT 

33 

PERIOD 

OCT 

73 

COMMA 

E-13 

a* 


CHARAC  HTR  ** 

D I MENS  HTR  ** 

ADOER  HTR  •• 

HlO  HTR  10 
DEC10  OEC  10*0 
HI  HTR  1 
PROD  HTR  ** 

AMASK  OCT  77777 
FLOAT  OCT  233000000000 
NUMB  HTR  *« 

XLOC  HTR  H 
CORE  OCT  100001 
JJ  SYN  XI 
END 


•  *  i 


OEG  COMPUTER  DATA  SUBMITTAL  FORM 


Submitted  by:  _  D«te:  _ 

Program  No.  _  Est.  Time  _  Classification 

Special  Instructions: _ 


Address 

Value 

Address 

Value 

Address 

Value 

Address 

Value 

mm 

— 

Bl 

B 

■ 

B 

B 

i 

mm 

1 

NOTES: 


1.  A  value  of  aero  must  be  entered  as  0,  not  left  blank. 

2.  Decimal  pta.  may  be  omitted  if  understood  to  follow  the  rightmost  digit. 

3.  The  value  3  X  10"*  may  be  entered  as  .00003  or  3-5.  not  as  3  X  10. 

4.  The  factor  portion  of  a  value  may  not  contain  more  than  •  digits. 

5.  The  exponent  portion  of  a  value  must  lie  within  the  ratye  ±39. 

6.  Bxpooeots  may  be  omitted  if  sero.  K  not,  they  must  be  signed. 

7.  Mask  cards  should  be  indicated  by:  - b - . 


E-15 

(REVERSE  BLANK) 
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